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Abstract 

Control of hyperfine interactions is a fundamental requirement for quantum computing architecture schemes 
based on shallow donors in silicon. However, at present, there is lacking an atomistic approach including critical 
effects of central-cell corrections and non-static screening of the donor potential capable of describing the hyperfine 
interaction in the presence of both strain and electric fields in realistically sized devices. We establish and apply a 
theoretical framework, based on atomistic tight-binding theory, to quantitatively determine the strain and electric 
^ held dependent hyperhne couplings of donors. Our method is scalable to millions of atoms, and yet captures the 

O 

Q strain effects with an accuracy level of DFT method. Excellent agreement with the available experimental data 
^ sets allow reliable investigation of the design space of multi-qubit architectures, based on both strain-only as well 
as hybrid (strain-bheld) control of qubits. The benehts of strain are uncovered by demonstrating that a hybrid 
^ control of qubits based on (001) compressive strain and in-plane (100 or 010) helds results in higher gate hdelities 

kO 

O and/or faster gate operations, for all of the four donor species considered (P, As, Sb, and Bi). The comparison 

between different donor species in strained environments further highlights the trends of hyperhne shifts, providing 

lO 

predictions where no experimental data exists. Whilst faster gate operations are realisable with in-plane helds for 
P, As, and Sb donors, only for the Bi donor, our calculations predict faster gate response in the presence of both 
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in-plane and out-of-plane helds, truly benehting from the proposed planar held control mechanism of the hyperhne 
interactions. 
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1. INTRODUCTION 


During the last few years, there has been signihcant progress HHa towards the realization of quantum 
computing architectures based on shallow donors in silicon (Si) |6l[7]. Several techniques have been explored 
to implement precise control of the nuclear or electronic spins through wave function engineering of donors 
using either electric helds [8] or strain helds [a HO]. At the core of such approaches, controlled manipulation 
of the donor hyperhne coupling is a critical component. Previous theoretical studies have been primarily 
focused on the electric-held-dependent Stark shift of the hyperhne interaction [B] for donors in Si, which 
is now well understood from both theory [TT1 - IT5] and experiments [IlH6l[I7|. In comparison, the strain- 
dependence of the hyperhne coupling is relatively less studied, despite ohering a promising alternative to 
manipulate the hyperhne coupling of donors. The presence of strain, in contrast to the use of an electric 
held, eliminates the possibility of ionization, as the control of the donor wave function is mechanical rather 
than electrostatic. Additionally strain can drastically reduce valley oscillations of exchange coupling [THl 
[TB] , which would play an important role in held control of qubits in strained environments. Recent progress 
towards atomically precise fabrication of donors in strained Si provides a testbed to demonstrate the 
advantages of strain in the realization of donor-based qubit devices pni- Whilst the previous studies 
have exclusively considered strain or electric held ehects on the quantum control of the donors, it is clear 
that through valley physics there is a subtle interplay between these two ehects. This work establishes 
a multi-scale theoretical approach to provide an understanding of the impact of strain and electric helds 
simultaneously present in the qubit devices, and predicts that such a hybrid quantum control scheme can 
open new avenues for architectures with faster single spin gates and spin-dependent tunneling read-out 
strategies. 

Existing theoretical studies of the impact of strain on the hyperhne interaction of donors have been 
based on either valley-repopulation model (VRM) derived from ehective-mass theory (EMT) [THl I2IH23] 
or density functional theory (DFT) [B]. While the VRM model was useful in providing a hrst-order 
description of the hyperhne shifts for small strain helds, it failed to explain the experimentally measured 
hyperhne reduction at large strain helds [9]. The DFT calculations for strained Si:P exhibited good 
agreement with the experimental measurements for an extended range of strain helds [9] , highlighting the 
importance of atomistic approaches. However, this method is limited to only few-atom systems, and is 
consequently unable to reproduce the donor binding energy spectra and provide a detailed picture of the 
wave functions [21|. Therefore the requirement for a theoretical framework with an atomistic accuracy 
accompanied by scalability to large-scale realistic systems remains a critical challenging problem. 

Our work hlls this theory gap by establishing a multi-scale atomistic tight-binding framework, which in 
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contrast to other approaches [IHl [2SH2Z| explicitly includes central-cell corrections and non-static dielectric 
screening of the donor potential. The multi-million-atom simulations for strain-dependent hyperhne inter¬ 
action are benchmarked to a high level against both ab-initio approaches and experiment. The electric 
held dependence of the hyperhne is accurately captured by demonstrating excellent agreement with the 
experimentally measured Stark shift data na for all of the four donor species considered (P, As, Sb, and 
Bi). A clear understanding of the inhuence of strain on the physical properties of the donor is presented 
in terms of underlying valley physics. The performance prospects of unstrained and strained Si substrates 
are explored, uncovering the benehts of strain for qubit devices, in particular by showing that a hybrid 
control of qubits based on ( 001 ) compressive strain and in-plane helds (100 or 010 ) results in higher gate 
hdelities and/or faster gate operations for all of the four donor species. Due to recent research interests 
for As, Sb, and Bi donors dSl im 128], we also present a comparison among diherent donor species in 
strained Si environments, further highlighting the trends of strain and electric held induced shifts in the 
hyperhne couplings, and providing predictions at large strain helds where no previous experimental or 
theoretical data exists. A novel scheme of two-dimensional hyperhne control in strained environments is 
explored based on electric helds from top and side gates. Whilst faster gate operations are realisable with 
the in-plane helds for P, As, and Sb donors, only for the Bi donor, faster gate response is predicted in the 
presence of both in-plane (100 or 010 ) and out-of-plane ( 001 ) helds, truly benehting from the proposed 
planar held control mechanism of hyperhne control. 


2. THEORETICAL FRAMEWORK 

In our TB approach, the Si bulk band structure is reproduced using a twenty-orbital (sp^d^s*) basis [29|. 
The donor atom is placed at the center of a large Si box (40x40x40 nm^) consisting of roughly 3.1 million 
atoms, and is represented by a Coulomb potential, U(r), which is screened by non-static dielectric function 
for Si and is given by: 


U(r) = -(1 + Aee-"" + (1 - A) ee"^" - e”^") (1) 

er 

where e is the electronic charge, and the previously published values of e. A, a, f3, and 7 are taken from the 
literature [131 EH] • Recently we have demonstrated the importance of the central-cell corrections and the 
non-static dielectric screening of the donor potential to accurately reproduce the experimental Stark shift 
for the Si:As donors |I3|. We now show that the non-static dielectric screening is also crucial to accurately 
reproduce the experimentally measured strained hyperhne interaction at large strain helds. Therefore this 
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work extends the TB model using the non-static screening function to P, Sb, and Bi donors. The donor 
potential is truncated to Uq at the donor site, r=0. The values of Uq are adjusted to reproduce the 
experimentally measured binding energy spectra of the donors [31]. By using the Uq values of 3.5 eV, 2.2 
eV, 3.8365 eV, and 4.6668 eV for P, As, Sb, and Bi donors respectively, we calculate the binding energies 
of the ground states Ai within 1 /reV and the binding energies of the excited states (T 2 and E) within 1 
meV of the experimental values for all of the four donor species. It should be emphasized that whereas 
multi-valley EMT theories have been successful in htting the ground state binding energies [T3|, the htting 
of all of the three Is states simultaneously with such a level of accuracy has been inaccessible. The accurate 
htting of the excited state energies, as achieved in this work, is critically important for strain-dependent 
hyperhne studies, where the excited state E mixes with the ground state Ai as a function of the strain. 


3. RESULTS AND DISCUSSION 

For the study of the strain-dependent hyperhne interaction of donors, we hrst benchmark our model 
against the recent experimental data set for the P donors [9]. A commonly adopted procedure to induce 
strain is by using the lattice mismatch technique, where two materials with diherent lattice constants are 
grown on top of each other. Such a technique is depicted in Figure [^a), where P doped Si is shown on 
top of a Sii_a;Gea; virtual substrate. The amount of strain in the Si region can be tuned by varying the Ge 
fraction [x) in the substrate. The lattice constant of Sii_ 3 ;Ge 3 ; is larger than Si, and therefore induces a 
tensile strain in P-doped Si region along the in-plane directions (a|| > asi). Gonsequently, the out-of-plane 
lattice constant (aj,) shrinks in accordance with the Poisson’s ratio, leading to a compressive strain along 
the growth direction (a_L < asi). For a (OOl)-oriented Sii_a;Gea:/Si system, the growth axis is the z-axis 
and the growth plane is the xy-plane, implying that the z-valleys (xy valleys) will primarily experience the 
effect of a compressive (tensile) strain. 

3.1 Characterising strain effects through valley physics: 

Figure [^b) schematically illustrates different effects on the splittings of Is donor states. The valley con- 
hguration of the lowest energy ground state is also included. A simple effective mass theory without 
multi-valley effects, such as presented by Kohn-Luttinger [32|, would lead to a six-fold degenerate Is state 
as shown in Figure [^bl). In reality, the effect of valley-orbit interactions results in a splitting of the Is 
energies into three sets (Figure [^b2)). The lowest ground state Ai is a singlet state, which is made up of 
all six valleys with a conhguration of 1, 1, 1, 1, !}• The hrst triply-degenerate excited state (T 2 ) has 
the following valley conhgurations: T 2 a;=;^{l, -1, 0, 0, 0, 0}, T 2 y=-^{ 0 , 0, 1, -1, 0, 0}, and T 2 z=;^{ 0 , 0, 
0, 0, 1, -1}. The second doubly-degenerate excited states (E) are composed of the following valley conhg- 
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(a) Schematic of P donor in strained Si (b) Schematic of Is energies, degeneracy levels, 
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(d) Fractional change in hyperfine interaction of Si:P 
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FIG. 1. Benchmarking tight-binding theory against experiment and DFT method for P donor: (a) 

An artist’s view of a single P donor in the strained Si. (b) Schematic diagram to indicate the impact of (b2) 
valley-orbital interaction and (b3) strain on the splitting of Is energy levels and the valley configuration of the 
lowest ground state Ai. The degeneracy of each energy level is also labeled, (c) The energies of the lowest few 
states of the P donor in the strained Si as a function of the valley strain (y), with (cl) for a large variation of strain 
and (c2) the low strain region (—4 < y < 0) is zoomed-in to highlight the valley repopulation regime. (c3) The 
energy difference of the lowest two states as a function of y. (d) The computed fractional change in the hyperfine 
interaction A(y)/A(0) is plotted as a function of y, and is compared with the previously published experimental 
measurements, and calculations based on valley-repopulation model and DFT. 

urations: Ea,y=^{l, 1, -1, -1, 0, 0} and 'Ejxyz=^^{-^, -1, -1, -1, 2, 2}. The influence of strain on the donor 
energies can be understood in terms of their valley configurations: the valleys in the direction of compres¬ 
sive strain experience a reduction in energy (higher population) and that valleys in the direction of tensile 
strain exhibit an increased energy (lower population) [21]. Since the excited states consist of assymetric 
valley contributions, they experience different effects of strain and therefore do not remain degenerate in 
the presence of strain (Figure [^b3)). In our case, the tensile strain along the x and y directions will push 
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the states with xy-valley conhgurations (T 2 a;, T 2 J,, '^xy) up and the compressive strain along the z-axis 
will shift the states with z-valley conhgurations (Ai^ and T 2 ^) downward on the energy scale. The valley 
repopluation effect for Ai^ has also been shown by illustrated by showing z-valleys (indicated by the red 
color) larger in size when compared with the xy valleys (indicated by the green color), in the schemetic 
diagram of Figure [^b3). 

The impact of strain on the donor states can be characterized either directly in terms of the Ge fraction 
X in the substrate, or can be described in terms of a dimensionless parameter so called the valley strain, 
X, which was derived by Wilson and Feher [^, and is given by: 


X = 


3A, 


USi ~ UCe 
USi 


2Ci2 

1 + - I X 


C 
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( 2 ) 


Here the value of the uniaxial strain parameter is 8.6 eV, Cn and C 12 are the elastic constants of Si 
and the value of their ratio C 12 /C 11 is 2.6, 6Ac=12.96 eV is the energy splitting of the singlet (Ai) and 
doublet (E) states for the unstrained bulk P donor, asi=0.5431 nm and aGe=0.5658 nm are the bulk Si and 
Ge lattice constants respectively, and x is the concentration of Ge in the virtual Sii_a;Gea: substrate. For 
this study, we vary x between 0 and 0.5, which corresponds to a variation of x from 0 to ~ —49. For each 
value of X, a strained TB Hamiltonian [33] is solved to compute the strained donor energies and states. 

Figure [^cl) plots the P donor energies calculated from TB simulations for a large variation of the 
valley strain x- The valley repopulation effect primarily occurs for small magnitudes of x, so Figure [^c2) 
presents zoomed-in version of the plot for —4 < x < 0 to highlight this effect. The effect of applied strain 
is on donor energies is accurately captured by the TB theory, indicating a partial lift of the degeneracy of 
the T 2 states, splitting them into a single T 2 ^ state whose energy decreases due to compressive strain along 
the z-axis, and a pair of degenerate T 2 a; and "l 2 y states with their energies increasing due to the effect of 
tensile strain. The T 2 X and T 2 y states remain degenerate as the same magnitude of strain is applied along 
the X and y-axis (aa;=ay=a||). The strain completely lifts the degeneracy of the E states. The energy of 
the Yjxy state increases due to the effect of tensile strain. 

The remaining two states Ai and 'Ejxyz with contributions from all of the six valleys experience the effect 
of both compressive strain along the z-axis and the tensile strain along the x and y-axis, and therefore 
exhibit a nonlinear dependence on x- The strain mixes the E^.^^ excited state into the ground Ai state. 
For —4 < X < 0, the Ai state experiences the competing effects of the tensile and compressive strains. 
Initially the increase in a|| is larger than the decrease in aj^, so the energy of the Ai state slightly increases. 
However, at the same time strain depopulates x and y valleys and increases z-valley contribution. This 
reverses the change in Ai due to the effect of decrease in a.± being increasingly dominant on the increasingly 
z-valley like Ai state. For x < —5, the lowest energy state is dominantly a z-valley state, with mixing from 
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FIG. 2. Theoretically predicted trends among P, Sb, As, and Bi donors : The plots of (a) the Ai energies 
and (b) the fractional change in the hyperfine couplings A(x)/A(0) are shown as a function of the Ge content x. 

the Exyz state to form a new ground state Ai^. The E^yz state is primarily composed of x and y valleys, 
and therefore its energy increases with strain. It is noted that the T 22 state does not mix with the Ai^ 
state as it is composed of two z-valleys with opposite signs. 

A recent study HI has dehned the ionization held as proportional to the energy splitting ((5E) of the 
ground state Ai and the higher excited state 2po, which is roughly 34.1 meV for P donors in bulk Si [31]. As 
evident from Figure[^cl), the strain reduces 5E, which becomes roughly 24 meV for y —20, the predicted 
strain held to suppress the valley oscillations of the J-coupling between P donor pairs [IHlIini- This implies 
a reduction in the ionization helds for the strained P donors, which would be useful for recently proposed 
spin-dependent tunneling read-out schemes [3l|. The energy diherence between the lowest two donor states 
is relevant in estimating time scales which determine the adiabatic condition in time-dependent processes 
driven by the gate potential variation. Figure [^c3) plots this energy diherence (AE 12 = T 22 —Ai^) as 
a function of y, indicating a reduction in its value due to the mixing of ground and excited states. At 
y —20, we calculate AE 12 ~ 2.5 meV which is smaller than the EMT value of —3.3 meV |18j . 

3.2 Hyperfine control by strain: 

The hyperhne interaction A(0) is directly related to the charge density at the donor site, |'0(O)p. Only the 
Ai state has a non-zero charge density at the donor site, thus only this state contributes in the determination 
of A(0). The excited states T 2 and E do not contribute to A(0). The applied strain reduces the hyperhne 
coupling due to the following reasons: (1) Valley Repopulation Effect: strain removes contributions from 
the X and y valleys, and increases z-valley contribution of the Ai state due to mixing of the E^y^ state. The 
E^yz state does not contribute in A(0), so A(y) becomes less than A(0). (2) Crystal Deformation: strain 
deforms the crystal and changes the bond-lengths from their bulk unstrained values. This modihes the 
radial distribution of the donor states, which are scattered over several Si lattice sites around the donor 
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atom, leading to a reduction in hyperfine. 

The VRM model based on the EMT theory only considers the hrst effect and the reduction in A(x) is 
represented by an analytical expression derived in Ref. m- 

A{x) 1 

R(0) 2 

Based on this model, the fractional change in the hyperfine A(x)/A(0) as a function of y is plotted in 
Figure [^d) using a dashed green line, along with the experimentally measured values (red dots) from 
Ref. [9]. Although the VRM method successfully describes A(y)/A(0) for small values of strain, it fails 
to capture the strain effects for the larger values of the applied strain (y < —10). In fact the VRM 
model limits the value of A(y)/A(0) to 1/3 for y < —20, based on the fact that all of the six valleys 
have equal contributions to the Ai state, whereas only two z-valleys contribute to the Ai^ state and hence 
A(y)/A(0)=2/6. Adding radial redistribution effects in the VRM model only reduces A(y)/A(0) to 0.3 
for y=—89 [22], still considerably different from the experimental data shown in Figure [^d), indicating 
A(y)/A(0) already reduced to ~ 0.22±0.09 at y ~ —29.5. Therefore in order to fully understand the strain 
dependence of the hyperfine, a more complete theoretical approach is required which takes into account 
both the valley-repopulation effect, as well as the volume deformation effect at atomistic scale. Recently 
reported DFT simulations [9] confirmed this notion by exhibiting a good match with the experimental data 
for both small and large values of strain (diamonds in Figure [^d)). Our TB calculations of A(y)/A(0) 
as a function of y are shown in the Figure [^d) using the square symbols, which demonstrate an excellent 
agreement with the experimental data as well as with the DFT calculations. For example, at y ~ —29.5, 
we calculate A(y)/A(0) as ?«0.284, compared to the experimental value of f»0.22±0.09 and the DFT value 
of ^0.27. It is noted that the previously applied static dielectric screening of the donor potential in the TB 
approach [T^ 125] 127] 135] results in a signihcantly higher value of «i0.364 for A(y)/A(0), which emphasizes 
on the requirement of the non-static (fc-dependent) screening of the donor potential for the study of the 
strain effects. The successful benchmarking of the TB method is in particular useful, because this approach 
has many advantages over the continuum EMT model and the computationally restricted DFT method. 
The TB theory not only accurately captures the atomistic physics, it is also scalable to simulation domains 
with millions of atoms, thereby enabling investigation of multi-qubit architectures. 

3.3 Tight-binding predictions for As, Sb, and Bi donors: 

After benchmarking the TB theory against the experimental data set for the Si;P system, we apply it to 
predict the influence of strain for other three donors (As, Sb, Bi), which have drawn signihcant recent 
research interests [161 HUES]. Figure a) plots the energies of the lowest donor state as a function of 
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the Ge concentration x in the Sii-a^Ge^; substrate, which is a more relevant parameter for experimentalists 
compared to x rnainly used in analytical theories. Overall the changes in the Ai^ energies follow similar 
trends for all four donors species. Figure [^b) plots the strain-dependent hyperhne A(a;)/A(0) as a function 
of the Ge fraction x. Again the overall trends are same for all the donors: in the valley-repopulation 
regime, the strain-dependent hyperhne decreases sharply but for larger strain values, it becomes much less 
dependent on the applied strain. Interestingly, for a given value of x, the order of AA(x)=A(x)/A(0) 
follows the same trend as the Stark shift parameter 'r] 2 - AA®’’(x) < AA^(x) < AA^®(x) < AA®‘(x) which 
is same order as One would naively expect this order to depend on the absolute 

values of the hyperhne interactions (A^(0) < A®’’(0) < A^®(0) < A®‘(0)), which is not true and in fact this 
sequence is directly related to the order of the binding energies of the donors (Af'’ < Af < Af^® < Af^). 
3.4 Benchmarking Stark shifts of hyperhne couplings against experiments: 

We have hitherto discussed the strain ehects on the hyperhne interactions of donors, which will be useful for 
the psoposed all-mechanical control of qubits uni. However, an alternative quantum computing architecture 
scheme could be based on a hybrid control of qubits, where the donors are present in strained Si and the 
control is applied by electric helds. Such hybrid control mechanism has certain benehts over traditional 
unstrained Si based systems, as the applied strain is expected to reduce valley oscillations of exchange 
coupling [UlIlH], as well as ionization helds. To study the ehects of strain and electic helds simultaneously 
present in the qubit devices, we hrst calculate and benchmark the Stark shift of the hyperhne interactions for 
all the four donor species under study. The Stark shift is calculated by adding a potential corresponding to 
an electric held of magnitude varying from 0 and 0.5 MV/m in the diagonal elements of the TB Hamiltonian, 
and computing the hyperhne interaction A{E) from the held-dependent ground state Ai. The change in 
the hyperhne relative to the absolute value of the hyperhne is then htted to a quadratic held dependence 
given by Eq. and the Stark shift parameter ( 772 ) is computed by htting to our simulation data [121 US]- 


A (b) - A (0) 
A(0) 


ri2\E 


(4) 


The computed values of r ]2 are plotted in Figure as a function of the Ai binding energies, which are in 
good agreement with the available experimental measurements. 

3.5 Hybrid (Strain-|-Field) control of qubits: 

For a bulk donor in unstrained Si, the x, y, and z axes are equivalent with respect to the application of 
an electric held. However for the donors in strained Si, where a|| 7 ^ a±, the ehect of an electric held along 
the in-plane direction is expected to be diherent from its ehect along the growth direction. Based on this 
notion, we suggest a new method of two-dimensional control of the donor hyperhne interaction by applying 
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FIG. 3. Comparison of the calculated Stark shifts with experiments: The calculated values of the Stark 
shift parameter rj 2 are plotted as a function of the corresponding ground state binding energies (Ai) of shallow 
donors in bulk Si, which are in good agreement with the available experimental values. 

electric fields along both the growth direction (z-axis) and one of the in-plane directions (x or y-axis). Such 
a scheme is schematically illustrated in Figure |^a), where the P-doped strained Si is displayed on top of 
Sii_a;Gea; substrate. The direction of electric held can be controlled by top and side gates, which apply 
helds E± and along the growth (001) and in-plane (010) directions, respectively. In atomically precise 
structures, in-plane side gates are frequently used to create an electric held {E\\) across the device |21 ES] 

. The impact of E\\ and Ej_ helds is investigated by varying their magnitudes from 0 to 0.5 MV/m, and 
calculating the Stark shift r] 2 {x) for each value of the strain, characterized in terms of x. 

Figure l^b) and (c) plots the strain-dependent Stark shift values {r] 2 {x)) as a function of x for all of the 
four donor species, when the applied held is (b) E\\ and (c) E±. In both cases, the magnitudes of ri 2 {x) 
increase overall as a function of the strain. The increase in \r] 2 {x)\ is larger for E\\ compared to Ej_ for the 
same x. This is because of the fact that the strain compresses the spatial distribution of the donor wave 
function along the z-direction and therefore the ehective Bohr radius along the z-axis is smaller than its 
values in the in-plane directions inn. Consequently for the same magnitude of the electric held, the net 
Stark ehect is stronger for the E\\ helds than for the E± helds. 

3.6 Strain leads to faster gate operations: 

In a quantum computer placed in a static held, donors of the same species will lie at or close to resonance. 
The ability to Stark tune an individual, targeted, spin away from this resonance allows for the addressability 
of an individual qubit. However, due to a hnite linewidth of an excited transition. Stark tuning by a larger 
frequency leads to higher hdelity and/or faster gates are achievable at the same hdelity. The timescale of 
an individual spin rotation is limited by the change in frequency provided by the Stark shift [6l [H]: Af{E, 
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(a) Schematic of planar control of 
hyperfine interactions 




FIG. 4. Planar control of hyperfine interactions in strained environments: (a) The schematic diagram 
illustrates a two dimensional control of the hyperfine interaction by applying E± and Fi|| fields along the z and 
y axes respectively. The plot of the strained hyperfine Stark shift r] 2 {x) for (b) the in-plane field E|| and (c) the 
out-of-plane held E_l as a function of the applied strain (x). 

x) = 'r] 2 {x)E‘^A{x)mi, where mi is the nuclear spin quantum number. We are interested in comparing the 
performance prospects of devices based on the strained Si with that of the unstrained Si, therefore we only 
compare the ratio of the two cases: Af{E, x)/Af{E, 0) = ri2{x)A{x)/ri2A{0). It is interesting to note that 
for the strained donors, the A(a;) value decreases (see Figure |^b)) as a function of the strain, but \ri 2 {x)\ 
increases (see Figures]^ (b) and (c)). A multiplication of these two quantities leads to Af{E, x)/Af{E, 
0) > 1 for E\\ helds, with its largest value being approximately 3 for the Bi donor. Therefore the E\\ held 
allows for faster control of the qubits in the strained Si substrate in comparison to the unstrained Si. On 
the other hand, Af{E, x)/Af{E, 0) < 1 for the Ej_ helds for the P, As, and Sb donors, thereby implying 
slower operation with the same hdelity for these donor species. Only for the Bi donor, we calculate that 
both the E\\ and E± helds exhibit Af{E, x)/Af{E, 0) > 1, which is promising given that the decoherence 
times for Si:Bi have been reported as comparable to that of the Si:P [38] . 


4. CONCLUSIONS 

A multi-scale atomistic framework, established by explicitly describing the cental-cell corrections and 
the non-static dielectric screening of the donor potential, was used to quantitatively study the individual 
as well as simultaneous ehects of strain and electric helds, as present in various quantum computing 
architecture schemes currently in development. Our calculations were based on millions of atoms in the 
simulation domain, and yet described the experimentally measured strain induced hyperhne shifts with an 
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accuracy level of the DFT method. We showed that a hybrid control scheme, where the donors are placed 
in a strained environment and the control of qubits is by electric helds, offers several advantages such as 
lowering of the ionization energies/helds and an increased magnitude of the Stark shift from the in-plane 
electric fields leading to the higher fidelity of single spin gates for all of the four donor species considered. 
The work demonstrates that the application of both strain and electric helds, and understanding their 
subtle interplay, has important implications for quantum control in the implementation of Si-dopant based 
quantum computing architectures. 
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